Experimental proposal for symmetric minimal two-qubit state tomography 
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We propose an experiment tiiat realizes a symmetric informationally complete (SIC) probability- 
operator measurement (POM) in the four-dimensional Hilbert space of a qubit pair. The qubit pair 
is carried by a single photon as a polarization qubit and a path qubit. The implementation of the 
SIC POM is accomplished with the means of linear optics. The experimental scheme exploits a 
new approach to SIC POMs that uses a two-step process: a measurement with full-rank outcomes, 
followed by a projective measurement on a basis that is chosen in accordance with the result of the 
first measurement. The basis of the first measurement and the four bases of the second measurements 
are pairwise unbiased — a hint at a possibly profound link between SIC POMs and mutually unbiased 
bases. 



Quantum state tomography, the procedure for infer- 
ring the state of a quantum system from measurements 
applied to it, is an important component in most, if not 
all, quantum computation and quantum communication 
tasks. The successful execution of such tasks hinges in 
part on the ability to assess with high efficiency the state 
of the system at various stages. 

A general measurement in quantum mechanics is a 
probability-operator measurement (POM). A POM is in- 
formationally complete (IC) if any state of the system is 
determined completely by the probabilities for the POM 
outcomes jll-[3]. State tomography infers these probabil- 
ities from the data acquired with the aid of the POM. 

A symmetric IC POM (SIC POM) is an IC POM of 
a particular kind. In a d-dimensional Hilbert space (of 
kets) it is composed of (P subnormalized rank-1 projec- 
tors, {'Pjl^^i, with equal pairwise fidelity of l/{d + 1). 
Their high symmetry and high tomographic efficiency 
have attracted the attention of many researchers, and 
a lot of work, both analytical and numerical, has been 
devoted to the construction of SIC POMs in various di- 
mensions, see e.g. [3-Q- 

In contrast to the major theoretical progress, up to 
date, all experiments and even proposals for experiments 
implementing SIC POMs have been limited to the very 
basic quantum system, the two-level system (qubit) [lOJ, 
ini, with the exception of ^U] where a SIC POM for a 
three-level system was approximated. This is, in part, 
due to the fact that there is no systematic procedure 
for implementing SIC POMs in higher dimensions, in a 
simple experimental set-up. 

In this contribution we suggest a feasible experiment 
that implements a SIC POM for four-dimensional Hilbert 
space of a qubit pair. The qubit pair is carried by a single 
photon, and the measurement is realized by passive linear 
optical elements and photodetectors. The experiment is 
clearly feasible with current technology. 

The proposal exploits a new theoretical approach for 
constructing SIC POMs JJj]. In this approach, we 'break' 
the SIC POM into two successive measurements, each 
with d outcomes, with the intention that each measure- 
ment would be relatively easy to implement. Unexpect- 



edly, we find that the two successive measurements that 
make up the SIC POM in dimension 4, not only have a 
strikingly simple form, but also hint at a close relation 
between the structure of mutually unbiased bases (MUB) 
and that of the SIC POM; Ref. [l^ is a recent review on 
MUB. 

Let us begin by setting up nomenclature and notations. 
A general measurement on a quantum system is com- 
posed of a set of outcomes. The latter are mathematically 
represented by positive operators Vj that sum up to the 
identity operator. The probability of obtaining the out- 
come Vj is given by the Born rule: pj = tr {Vjp), where p 
is the pre- measurement statistical operator of the system. 
If the jth outcome is found, the post-measurement sta- 
tistical operator of the system is given by pj = —PjpP^, 
where Pj is the relevant Kraus operator for the jth out- 



come, Vi 



P'P 

^j n ■ 



Note that the decomposition of the 



Vs into the corresponding Kraus operators is not unique; 
for example, P]Pj is invariant under the unitary trans- 



formation Pj — >■ C/,Pi 



jjj, with different C/jS corresponding 



to different implementations of the POM. 

Suppose that a given system is subjected to a sequence 
of two POMs, each with d outcomes, {Ak — ^I^fclfe^ii 
followed by {Sj'°'}^^i, where the superscript k indicates 
that in general the second measurement depends on the 
actual outcome of the first measurement. Following 
Born's rule, the probability of obtaining the nth and mth 
outcomes for the first and second measurements is given 
by tr (^pA^^B!^^ An) ■ Accordingly, the two successive mea- 
surements are equivalent to a single POM with d^ out- 
comes Vn.m — Al^B'^'An with 71,771—1, . . . , d. Indeed, 
upon finding the over-all outcome Vn,rm we know that 
the nth outcome of the first POM and the mth outcome 
of the second POM are the case. In what follows we will 
identify the ^s and the Bs such that the Vs make up 
the SIC POM for a qubit pair. But before doing so, we 
discuss such an identification for a single qubit. 

All SIC POMs in two-dimensional Hilbert space are 
unitarily equivalent to the "tetrahedron measurement" 
(TM), whose outcomes correspond to four vectors that 
define a tetrahedron in the Bloch sphere 0, [^] . The TM 
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FIG. 1: An optical implementation of the tetrahedron mea- 
surement using two successive measurements. 



could be realized by a sequence of two measurements, as 
sketched in Fig.[T] Here, the qubit is encoded in a spatial 
alternative of a single photon ("path qubit"): traveling 
on the left or on the right. A unitary transformation on 
the qubit state amounts to sending the photon through a 
set of beam splitters (BSs) and phase shifters (PSs) [l^ . 
In this optical setting the TM is implemented as 
follows: First, two BSs (BSl and BS2) are used to 
implement the Kraus operators v4i = diag(ii, ^2) and 
A2 = diag(ri,r2), where diag stands for a diagonal ma- 
trix, and ti and r^ are the transmission and reflection 
amplitudes of the ith BS. A photon which enters the 
apparatus with a path statistical operator p, exits at 

port k with the statistical operator AkpAytr (AkpAj. 



For the values ty- 



these operators correspond to the measurement outcomes 
^^ = 1(1 -I- 1 (_i)fecr3) with fc=l, 2 (the cts are the Pauh 
operators with (T3 diagonal in the left-right basis). Then 
a photon that exits the first measurement apparatus at 
port 1 is measured in the ai basis while a photon that 
exits at port 2 is a measured in the cr2 basis. These mea- 
surements could be realized by balanced BSs and appro- 
priate PSs, as indicated in the figure. 

Actually, the TM was successfully implemented in an 
optical system [I^l, where the qubit was encoded in a 
photon's polarization ( "polarization qubit" ) rather than 
in spatial alternatives. The set-up of [ifll also consisted of 
a sequence of two measurements, quite analogous to what 
is described above. In that set-up, a partially polarizing 
beam splitter (PPBS) was used to implement the Kraus 
operators Ak and then, depending on whether the photon 
was transmitted or reflected, a measurement of cti or (72 
followed. 

Before moving on to the qubit-pair case, let us close the 
present discussion with three remarks: (i) In the above 
construction, the qubit MUB play a central role; they are 
used to construct, by means of successive measurements, 
the SIC POM. We will see below that such a relation 
appears in dimension 4 as well, (ii) A practical imple- 
mentation of the scheme presented in Fig. [1] requires the 
stabilization of the interferometer loop defined by the 
four BSs. (iii) SIC POMs for a three-level system could 
be implemented by using a similar set-up, but with al- 



lowing the photon to take three different paths [13]. 

In dimension 4, there is only one known SIC POM, and 
all the other known SIC POMs are unitarily equivalent 
to it 6]. This SIC POM is composed of 16 subnormal- 
ized projectors onto 16 (fiducial) kets. The latter are 
represented in the following matrices as columns with 

N = v/sTTH and x = V^ + V^ 111, 




Each of these matrices could be written as a diagonal 
matrix times a unitary matrix. The set of bases, cor- 
responding to each unitary matrix, together with the 
computational basis, form the complete set of MUB in 
dimension 4. To be more specific, the diagonal matrices 
are 

Ai =—diag(x, 1,1,1), A3 = — diag(l,l,x, 1), 

A2=^diag(l,x,l,l), A4 = ^diag(l,l,l,x), (2) 

and the unitary matrices are 
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(3) 



Noting that ^ ■ A!:Aj = 1, we identify the ^s with the 
Kraus operators of a measurement. 

Actually, the operations of Eq. (|3]) transform the com- 
putational basis into the MUB, 
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FIG. 2; A successive-measurement scheme for realizing the SIC POM of a qubit pair. Here the two-qubit state is encoded in 
the spatial-polarization state of a single photon. 



where CZ stands for the controUed-Z (phase flip) oper- 
ation, CZ — diag(l, 1, 1, —1). The bases fBi and ^82 are 
composed of product states, while the bases QSa and ©4 
consist of maximally entangled states. 



The structure of fiducial vectors in Eq. ^ (which form 
the SIC POM in dimension 4) allows us to implement the 
SIC POM by two successive measurements: A measure- 
ment whose Kraus operators are given in Eq. ([2]), and 
depending on the measurement outcome, a measurement 
in one of the MUB of Eq. Q. Next, we propose an op- 
tical implementation for this scheme. 



Our proposal is based on the methods of [15'] where 
the two qubits, a polarization qubit and a path qubit, 
are encoded in a single photon. (We chose here to use 
a polarization qubit instead of another path qubit in or- 
der to avoid as many interferometric loops as possible in 
the optical set-up.) We consider the vertical (v) and hor- 
izontal (h) polarizations as the basic alternative of the 
polarization qubit, and traveling on the left (L) or on the 
right (R) as the basic alternative of the path qubit. A 
unitary transformation on the two-qubit state amounts 
to sending the photon through a set of passive linear op- 
tical elements that unitarily change the state of the path 
and polarization qubits |15| . For our purpose we need 
the following optical elements: half- wave plates (HWPs), 
BSs, polarization and path dependent PSs, and PPBSs. 
A PPBS is a BS whose reflection and transmission coeffi- 
cients depend on the polarization. Its action corresponds 
to a joint unitary transformation on the polarization- 
path qubits. In the present context, it suffices to con- 
sider a PPBS with real amplitude division coefficients r 
and t that obey the unitarity condition r^ +t^ = 1 for 



the vertical and horizontal polarizations. 
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(5) 



This is a block-diagonal matrix, with the blocks trans- 
forming the vertical or horizontal polarization, respec- 
tively. Two cases of interest are (i) r^, = th = !'■ the po- 
larizing beam splitter (PBS) which totally reflects (trans- 
mits) vertically (horizontally) polarized hght, and (ii) 
fv —Th — 1: the CZ gate. 

The Kraus operators for the first measurement are 
listed in Eq. ^. Their realization is schematically 
drawn in Fig. [2] at the 'first measurement' part. For 
each port, we set the parameters of the different opti- 
cal elements such that a photon which enters the ap- 
paratus with a polarization-path statistical operator p, 
exits at port k with the two-qubit statistical operator 

AkpA^^./iY{AkpA\\. To be more specific, the appara- 
tus is configured such that the beam splitters BSla and 
BSlb have the same properties and so have beam split- 
ters BS2a and BS2b. The PPBSs on the left and right 
arms also have the same properties. The reflection co- 
efficient of BSla and BSlb is ri = 1/A^. The refiection 
coefficient r2 of BS2a and BS2b satisfies tir2 = 1/iV, that 
is, r2 — 1/VN'^ — 1, where ii is the transmission coeffi- 
cient of BSla(b). Setting r^ — th = y in Eq. ([SJ, the 
two PPBSs transform vertically polarized incident light 
|v) to the polarizations y |v) and \/l — y^ |v) in the re- 
fiected and transmitted arms, respectively, and horizon- 
tally polarized light |h) to the polarizations \Jl — y^ |h) 
in reflection and y |h) in transmission. The amplitude di- 
vision coefficient y is chosen such that tit2y = ^/N, and 
therefore, y = 1/^/N^ 



2, where ^2 is the transmission 



coefBcient of BS2a(b). These settings ensure that the 
measurement of Eq. ^ is reahzed. 

To complete the measurement scheme, a second mea- 
surement is taking place. This measurement depends on 
the actual outcome of the first measurement, namely, on 
the output port where the photon exits. For photons 
emerging from the fcth port, basis Sfc of Eq. ^ is mea- 
sured. In order to measure in a given basis, Sfc, we 
first apply a unitary operation Uk of Eq. ([3]) that trans- 
forms the basis 25^ into the computational basis and then 
measure in the computational basis by using PBSs and 
photodetectors, as illustrated in Fig. [2] at the 'second 
measurement' part. 

To implement the unitary transformations of Eq. ([3|). 
one could use either a single, specially designed, bircfrin- 
gent material, or a sequence of wave plates and PPBSs. 
Considering the latter option, these unitary transforma- 
tions are 

U2 = (C/psC/hwp) ® (UpsUus) , 
Z^3 = Cz((C/psC/hwp)'»C/bs), 



Ui^CzfUn^ 



(C/psC/bs) , 



(6) 



where C^s — diag(l,i) shifts the phase of the path and 
polarization qubits by n/2, and L^hwp and Ubs implement 
the Hadamard gate. 
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(7) 



for the polarization qubit and the path qubit, respec- 
tively. For this aim, we use a HWP with its major axis 
at an angle tt/8 to the optical axis, and a balanced BS. 
We see that the unitary transformation Ui and U2 
can be decomposed into a tensor product of two uni- 



tary transformations, one for the path qubit and one for 
the polarization qubit. The unitary transformation U3 
and W4 are not of that kind and could be realized, for 
example, by using PPBSs together with a Mach-Zehnder 
interferometer. This closes our proposal. 

To conclude, we are proposing a feasible experimen- 
tal scheme that implements the SIC POM for a two- 
qubit system. Our scheme uses linear optical elements 
and photodetectors, and is, therefore, well within the 
reach of existing technology. The proposal is based on 
a successive-measurement approach to SIC POMs. We 
found that the SIC POM for the qubit pair corresponds 
to a POM diagonal in the computational basis, followed 
by projections onto bases which are mutually unbiased. 
We observed that this unique construction is owed to a 
structural relation between the fiducial vectors and the 
MUB vectors in dimension 4. 

On a more general note, we believe that it would be 
interesting to learn, if and how this scheme can be gen- 
eralized to higher dimensions. Such a study could be of 
a theoretical and a practical use; it might teach us about 
the SIC POMs' structure in high dimensions and pro- 
vide new ideas for implementing them. For example, it 
is possible to show [l^ that any SIC POM which is co- 
variant with respect to the Heisenberg-Weyl group can 
be realized by two successive measurements, each of a 
rather simple form. This is, in particular, interesting for 
dimension three, where one has a one-parameter family 
of non-equivalent group-covariant SIC POMs [6]. 
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